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SUMMARY 


Electron concentration and collision frequency were evaluated at several 
axial flow stations along a blunt-nosed 9° half-angle cone flying at zero angle 
of attack, 17/700 ft/sec, and 170,000-foot altitude. These conditions were chosen 
to correspond to those for the maximum observed attenuation during the flight of 
the RAM A1 vehicle. Equilibrium, frozen, and finite-rate flow assumptions were 
made about both chemical and ionic processes. Boundary- layer concepts were applied 
to modify the results near the body. Signal attenuation at 240 me was computed by 
numerical integration through the plasma layer of the propagation equation for 
electromagnetic radiation. The propagation model was a monochromatic plane wave 
normally incident on a plane-parallel plasma in which electron concentration and 
collision frequency varied only in the propagation direction. 


INTRODUCTION 


The ionized layer around a space vehicle during reentry strongly attenuates 
radio signals and therefore often breaks communication links with the ground. The 
importance of the problem has been recognized for some time (ref. l) and many 
papers have been written which pertain to it in one way or another. 

Most of the experimental flight data on signal attenuation have come from 
telemeter records of flights which had objectives other than the study of communi- 
cations problems . As a result, these data commonly show little more than thqt 
the signal blacked out during a certain part of the trajectory. The RAM (Radio 
Attenuation Measurement) project at the Langley Research Center has obtained 
better data than those by flying blunt-nosed slender bodies designed to give meas- 
urable attenuation without signal blackout and to give quantitative information 
about the amount of signal attenuation, the effect of the plasma on the antenna, 
the plasma properties, and the attitude, altitude, and velocity of the vehicle. 
Reference 2 reports the results of the RAM A1 flight. 

A good theoretical understanding of the signal attenuation problem is neces- 
sary not only for intelligent interpretation of flight data but also for advance- 
ment of the art of designing future communications systems. There are two basic 
parts to the theoretical problem: the specification of the plasma parameters, and 

the calculation of the effect of the plasma on the electromagnetic wave which 
passes through it. 


For a hypersonic vehicle the first part of the problem involves a calculation 
of the real-gas flow field around the body in sufficient detail to give the elec- 
tron concentration and collision frequency profiles through the shock layer at the 
body stations corresponding to antenna locations. The profiles obtained depend 
strongly upon whether chemical and ionic reaction rates are fast or slow relative 
to flow time over the body. Important differences in the profiles are also pro- 
duced by allowing for a boundary layer near the surface. No calculations for com- 
plete nonequilibrium flow fields including boundary layer have appeared in the 
literature. However, several useful calculations for inviscid, real-gas, blunt- 
body flow fields have been reported. Typical papers from this group are refer- 
ence 3 , which reports thermodynamic and flow properties for equilibrium air in 
inviscid flow over a blunted 9° half -angle cone, and references 4, 5, and 6 , which 
report the calculated finite-rate variation of electron concentration along stream- 
lines in inviscid flow over some other blunt bodies. Calculations for chemically 
reacting boundary layers have been reported in references 7 , 8 , and 9 * 

In order to test the effect on calculated radio signal attenuation of various 
assumptions about the nature of the flow without a large expenditure of effort on 
computations, new methods have been developed, and are described herein, by means 
of which real-gas, nonequilibrium flow fields, including boundary- layer effects, 
can be specified with reasonable confidence from normalized shock- layer profiles 
derived from published results of detailed equilibrium inviscid calculations over 
a wide range of flight conditions. 

The second basic part of the theoretical problem requires the solution of 
the propagation equation for electromagnetic waves through a region in which the 
conductivity of the medium is complex, variable, and characterized by gradients 
which are too steep to permit WKB solutions (ref. 10) to be used. Analytic solu- 
tions valid for steep electron concentration gradients are available for special 
profiles (ref. 1 1 ), but these do not offer the flexibility necessary for general 
use. A unique method has been devised for accomplishing numerical integration of 
the propagation equations. Although numerical integration of these equations 
through the plasma in the direction of propagation is not possible because of the 
contribution of reflection to the initial boundary conditions, the equations can 
be integrated in the reverse direction if the transmitted electric vector is nor- 
malized to unit amplitude and zero phase angle. Solutions of this kind are used 
herein, and the method has been checked by comparison of special calculations to 
published results for profiles which lead to analytic solutions (ref. 12). Sub- 
sequent to the development and first publication (ref. 13 ) of the foregoing method 
(for normal incidence only), a paper (ref. l4) appeared which described the inde- 
pendent development of essentially the same approach to the problem of solving 
the propagation equations for arbitrary angles of incidence. A computer program 
for arbitrary angles of incidence was already in operation at the Langley Research 
Center when reference 14 came to the authors' attention, but the results were not 
published until some time later (ref. 12 ). 

A brief preliminary presentation of results was given at the Second Plasma 
Sheath Symposium (ref. 13) • Reference 13 also includes a comparison of calculated 
and observed attenuation during the flight of the RAM A1 vehicle. 
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SYMBOLS 


Cp 

D n 

E 

Er 

e t 

E = 
E R = 


pressure coefficient. 


P ~ Po 


nose diameter, 5.08 cm for RAM A vehicle 
electric field intensity, volts/meter 
electric field intensity reflected by plasma layer 
transmitted electric field intensity 
E 


Erp |e 


i ( k o A+ 0b) 


E- 


'R 


Erp I 0 


i(koA+^j 


electron 

real part of complex number representing E 


f = 


g 


M 


imaginary part of complex number representing E 


|E T 


-i( k o A+ 0b) 


I = 


N 

N2 

WO 

W 0 + 


i T 


i =fl 


wave number; reaction rate (dimension of wave number is cm~ E 
sion of reaction rate is cm^-molecule-E-sec -1 ) 


atomic nitrogen 
molecular nitrogen 

nitric oxide 
nitric oxide ion 


; dimen- 
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atomic oxygen 
molecular oxygen 
pressure, atm 

velocity-averaged effective collision cross section for electrons with 
neutrals in air, cm^ 

dynamic pressure, atm 

universal gas constant 

radius of body at x/D n , cm 

entropy 

signal loss due to reflection, db 
total signal loss, db 

distance along streamline from shock front, cm 
temperature, °K 
time, sec 

velocity, ft-sec~l 
real part of square of k/k D 
imaginary part of square of k/k D 
mass-flow parameter, slug-sec - -*- 

distance along body axis from stagnation point, cm 

distance from body along normal to body surface, cm 

degree of dissociation of air; phase parameter in equations (9) 
and (12), cirT^ 

attenuation parameter in equations (9) and (13), cm'-*- 
ratio of specific heats 

thickness of shock layer measured along y, cm 
shock standoff distance at nose, cm 
boundary- layer thickness, cm 
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collision frequency of an electron with neutral particles in air, 
density, slugs-ft - ^ 

phase angle of E relative to phase of EQ a , radians 
angular signal frequency, radians-sec~l 
angular plasma frequency, radians-sec - ^ 
concentration, particle cm - 5 
magnitude of a complex number 
Subscripts: 

0 medium containing incident and reflected electromagnetic waves 

1 plasma portion of wave medium 

2 medium containing transmitted electromagnetic waves 

a boundary between medium 0 and medium 1 

b boundary between medium 1 and medium 2 

bZe boundary- layer edge 

bs body streamline (in inviscid flow) 

d dissociation rate of molecules in air 

ns normal shock 

o free-space value 

r recombination rate of electrons in air 

s shock property 

sZ standard sea-level conditions 

oo free stream 

Flow-field notation: 

EQ flow field in which chemical and ionic reactions are everywhere in 

equilibrium 

EQBL same as EQ but with boundary- layer corrections to temperature and 

density 


0 


CUp 

n 

i i 


sec 


-1 
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FD flow field in which chemical reaction rates vanish downstream of freeze 

point on each streamline 

FDFKBL same as FD but with boundary-layer corrections to temperature, density, 
and velocity and also with allowance for finite rate recombination of 
electrons along streamlines 

FS flow field in which chemical and ionic reaction rates vanish downstream 

of shock wave, but shock is at equilibrium 

SFD same as FD except that £e~J is in local equilibrium 

SFS same as FS except that i s in local equilibrium 

Prime to a symbol denotes a derivative with respect to y. 

FLOW-FIELD CONSIDERATIONS 


A schematic diagram of the flow field over the RAM body in flight is given 
in figure 1. The coordinate system shown is used throughout the paper. 



Figure 1.- Flow- field diagram for RAM A vehicle. 

The equilibrium inviscid flow-field properties were determined by means of 
an interpolation technique which adapted numerical flow-field results (ref. 3 ) 
for equilibrium inviscid flow over a blunted 9° half-angle cone at selected 
velocity- altitude points to the velocity and altitude corresponding to peak 
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VHF signal attenuation during the flight of RAM Al. This technique is “based on 
the fact that entropy and pressure distributions across the shock layer can be 
normalized in such a manner as to be insensitive to real-gas effects (such as 
dissociation) over wide ranges of hypersonic flight conditions. (Parameters such 
as temperature and density are more sensitive.) The technique is similar to that 
used in classical boundary- layer theory in that the flow parameters are normalized 
in terms of their values at the boundaries of the shock layer (body and shock 
values) . 

The selection of entropy as one of the normalized parameters is advantageous 
because streamlines are isentropic lines in equilibrium inviscid flow. The selec- 
tion of pressure as one of the normalized parameters is conventional since shock 
and body pressure distributions obtained from Newtonian and blast wave approaches 
are usually expressed or compared to distributions from characteristics calcula- 
tions in this way. 

The normalized entropy and pressure distributions used are shown in figures 2 
and 3> respectively. Each curve corresponds to a different body station. The 
distributions represent the results of reducing numerical solutions for a wide 
range of real-gas flight conditions as, for example, those published in refer- 
ence 3 . The variation with flight conditions from the curves of figures 2 and 3 
is only a few percent. Distributions in reasonable agreement with those shown 
have also been obtained from other unpublished solutions. In fact, 7 = 1.4 
solutions at a Mach number of 8 also give distributions which are within the 
limits mentioned. 


X 



s ne " 


Figure 2.- Entropy distribution in equilib- 
rium inviscid flow over a blunted 9 ° 
half-angle cone. 



Figure 3 .- Pressure distribution in equilibrium 
inviscid flow over a blunted 9 ° half -angle 
cone. 
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In order to use the plots shown in figures 2 and 3, the body and shock values 
of entropy and pressure under the flight conditions of interest must be specified. 
It is in the insertion of these boundary values, particularly the body value, that 
the real-gas effects are reintroduced. This is done as follows: (a) entropy is 

constant along the body streamline and is equal to the normal shock or stagnation- 
point value at the nose; (b) entropy along the shock is related to normal shock 
entropy and to free-stream entropy by means of the normalized plot shown in fig- 
ure k; (c) body and shock pressures are related to the normal shock and free- 
stream pressures by means of the normalized pressure-distribution curves of 
figure 5. 


Figures 4, 5> and 6, like 
figures 2 and 3> represent the 
results of reducing real-gas 
numerical solutions such as those 
of reference 3. These plots are 
not as independent of flight con- 
ditions as are those of figures 2 
and 3" Thus, flight conditions 
most closely approximating those 
required were chosen for the prep- 
aration of these plots. Fig- 
ure 6 is a plot of shock-layer 
thickness along the body in terms 
of the shock standoff distance at 
the nose. When this curve is 
used, the value of the shock 
standoff distance can be found 
from the Lighthill equation: 





+ 0.002 



The real-gas flow properties 
behind the normal shock and for 
the ambient free stream as func- 
tions of altitude and velocity 
were taken from reference 15. 


Figure 4.- Entropy distribution along shock in 
equilibrium inviscid flow over a blunted 9° 
half -angle cone. 


The interpolation technique which has been described is useful as an expedi- 
ent means for approximating flow- field properties. It permits reasonably accurate 
evaluation of equilibrium inviscid flow-field properties over bodies of a given 
shape for arbitrary values of velocity and altitude over a fairly wide range. It 
is worth noting that the body values of the flow parameters can be more accurately 
determined than the shock values and that the flow parameters near the body are 
those of prime importance to radio-signal-attenuation predictions. The flow prop- 
erties computed for RAM A1 at the peak attenuation point on its trajectory are 
summarized in table I. The properties specified are pressure distribution and 
shock and streamline coordinates. 
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Figure 5*~ Pressure distribution along streamlines in equilibrium inviscid flow over a blunted 9° 

half -angle cone. 
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Figure 6.- Shock- layer thickness for inviscid equilib- 
rium flow over a blunted 9° half-angle cone. 


Shock coordinates, streamline coordinates, and pressure distributions were 
used without change for nonequilibrium inviscid as well as equilibrium inviscid 
flow, since these properties are relatively insensitive to changes in chemical 
composition for the RAM shape. Temperature, density, and composition for equi- 
librium inviscid flow were followed along streamlines by using pressure and 
entropy as inputs to real-air property tabulations. Temperature, velocity, and 
density for frozen inviscid flow were obtained from isentropic flow relations 
by using pressure distributions along streamlines and initial conditions at the 
freeze point as inputs. The specific heats and enthalpies used in the isentropic 
frozen-flow expansions were computed on the basis of fixed composition as deter- 
mined at the freeze point. The internal energies of the species were assumed to 
be in equilibrium with the translational temperature throughout the frozen 
expansion. 

Streamlines entering the boundary layer were followed by using a boundary- 
layer mass-flow parameter. The boundary-layer thermal properties and velocity 
distributions were determined by means of the laminar similar solutions treat- 
ment for equilibrium boundary layers described in reference 16. The way in which 
this boundary-layer treatment was adapted for use in frozen flow and also the 
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TABLE I.- ENTROPY AND PRESSURE ALONG STREAMLINES IN EQUILIBRIUM INVISCID FLOW OVER A 


BLUNTED 9 ° HALF-ANGLE CONE 
[Altitude = 170,000 ft; Velocity = 17,700 ft/secj 


(a) Entropy 

= 0.05; % = 51. 00; = 44.10 

R R 


Streamline 

s/r 

Body 

44.10 

0 

43.46 

1 

42.48 

2 

40.90 

3 

37.72 


(b) Pressure distribution and streamline coordinates 




Body 


0 


1 


2 


3 

Shock 

D n 
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( C p)bs 
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C P 

A 

( C P) S 
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( C P)ns 
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( C p)ns 
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( C p)ns 

A 

( C p)ns 
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( C p)ns 

Dn 

( C p)ns 

0.1 

0 

O.58O 

0.355 

0.647 

0.835 

0.750 





0.053 

0.770 

.2 

0 

.315 

.213 

.354 

.496 

.427 

0.915 

O.56I 



.079 

.588 

• 3 

0 

.167 

.168 

.187 

.418 

.241 

.660 

.313 



.114 

.450 

.42 

0 

.075 

.148 

.084 

.327 

.116 

• 579 

.164 

0.991 

0.320 

.172 

.325 

• 5 

0 

.064 

.147 

.069 

.317 

.086 

• 533 

.118 

.908 

.217 

.207 

.263 

•75 

0 

.046 

.134 

.047 

.288 

.052 

.500 

.065 

.807 

• 099 

.290 

• 153 

1.0 

0 

.038 

.119 

.038 

.259 

.040 

.446 

.046 

• 759 

. 0 64 

.349 

.110 

2.0 

0 

.026 

.082 

.026 

.183 

.026 

.329 

.027 

•593 

.031 

.504 

.065 

3.0 

0 

.022 

.059 

.022 

.140 

.022 

.266 

.022 

.501 

.024 

.601 

.049 

4.0 

0 

.020 

.050 

.020 

.119 

.020 

.231 

.020 

.446 

.021 

.663 

.041 

5-0 

0 

.019 

.041 

.019 

.103 

.019 

.206 

.019 

.405 

.019 

.705 

.036 

7-5 

0 

.018 

.034 

.018 

.082 

.018 

.165 

.018 

.340 

.018 

.774 

.030 

10.0 

0 

.019 

.029 

.019 

.073 

.019 

.146 

.019 

• 303 

.019 

.820 

.028 


method for matching the boundary- layer profile to the shock- layer profile are 
described later. 

As previously stated, the main objective of the study was to compare the 
effect on predicted radio signal attenuation of several assumptions about the 
nature of the flow. Three kinds of assumptions were made. First, the thermal 
state of the air was assumed to be completely determined by the neutral molecules 
and atoms. For the flight conditions chosen this is a valid approximation for the 
determination of thermodynamic properties, since the ion concentrations are too 
small to have an appreciable effect on the thermal properties. Calculations were 
made for the following limiting conditions: (l) complete chemical equilibrium 

throughout the flow, and (2) frozen chemistry subsequent to the attainment of 
complete chemical equilibrium in the shock front. A more realistic calculation 
was made in which dissociation occurred at a finite rate behind the shock front 
until a "freeze point" was reached. Typical streamline temperature distributions 
for these three conditions are shown in figure 7. 
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Figure 7 .- Temperature distribution on streamline 1 for equilibrium, frozen at shock, and frozen 

downstream cases. 


Second, assumptions were made about the ionic state of the air. The limiting 
conditions were: (l) .complete ionic (as well as chemical) equilibrium throughout 

the flow, and (2) frozen ions (as well as chemistry) subsequent to the attainment 
of complete equilibrium in the shock front. Two intermediate calculations were 
made on the basis of frozen flow downstream of the freeze point. In one of these 

the reaction N + 0 N0 + + e - was considered to be in equilibrium at and down- 

stream of the freeze point. In the other this reaction was assumed to be in 
equilibrium at the freeze point but downstream of this point a finite rate was 
used for the recombination of N0 + and e _ . 

Third, inviscid flow was assumed throughout the shock layer as a first 
approximation. Boundary-layer corrections were made only for the calculations in 
which complete equilibrium was assumed and for those in which finite-rate recom- 
bination of ions was assumed. Figure 8 shows the variation of electron concen- 
tration along a typical streamline for the following assumptions: 
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Figure 8.- Electron concentration along stream- 
line 1 for several chemical-ionic models. 

EQ equilibrium, complete thermal equilibrium throughout the flow 

FS frozen shock, complete freezing of composition at shock-front equilibrium 
values 

SFS semifrozen shock, same as FS but with N + 0 N0 + + e" in local 

equilibrium 

FD frozen downstream, same as FS except that freezing takes place downstream 
of shock and is not an equilibrium composition at the freeze point 

SFD semifrozen downstream, same as FD but with N + 0 ^ N0 + + e~ in local 
equilibrium 

FDFKBL frozen downstream with finite rate recombination of ions and with 

boundary- layer corrections, same as FD but with allowance for (l) a 
finite rate for N0 + + e“ -» N + 0 and (2) altered values of velocity, 
temperature, and density after the streamline enters the boundary layer 
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NONEQUILIBRIUM CONCEPTS 


In general the true plasma properties are best represented by calculations 
which include the effects of finite reaction rates. The most accurate way to do 
this is to use a computer program which simultaneously calculates composition, 
thermodynamic properties, and flow-field properties. The following approximate 
treatment of nonequilibrium processes was satisfactory for the purposes of this 
study . 

In this paper air is assumed to be a mixture of 21 percent oxygen and 79 per- 
cent nitrogen. The dissociation rate of initially undissociated air is controlled 
by the dissociation rate of Op, as can be seen from the calculations reported in 

references 17 and 18. A fairly accurate model for calculating the dissociation 
rate of air is therefore a dissociating diatomic gas for which the dissociation 
rate is identical to that of pure Op* However, because of the large difference 
between the dissociation energies of Op and Np, the thermodynamic properties 
of this gas are not accurately represented by dissociating Op. These thermo- 
dynamic properties were assumed to be identical to the properties of a dissoci- 
ating mixture which is initially 21 percent Op and 79 percent N2 and in which 
the ratio of 0 atom concentration to N atom concentration is fixed at 5 to 1 
and rotation and vibration are in equilibrium. The choice of this particular 

ratio was based on figure 9 which shows [®]/W as a ^ unc ti on J^OpJ behind a 

normal shock wave in air at shock temperatures comparable to the range of interest 
in this paper. (The curves in figure 9 are crossplots of figures 9 and 10 in ref- 
erence 17.) The value of [V]/[nJ for the 10,000° K shock temperature approaches 
5 to 1, but the 7, 000° K data call for a value of of the order of 50 

to 1. This factor-of-10 discrepancy is not as serious' as it seems at first sight, 
since the degree of dissociation produced along the cooler streamlines is small. 

The NO concentration was arbitrarily taken to be zero. This is equivalent 
to assuming that half the NO molecules are Op and the other half are Np. 

Since the specific heats of Op, NO, and Np are approximately the same, the 
most serious error involved is the neglect of the heat of formation of NO from 
O2 and Np, which causes the temperature predicted by the model to be too low. 

For example, if the mole fraction of NO is 0.1 and a = 0.2, then the tempera- 
ture error due to this cause is of the order of 300° K. 

The variation of the degree of dissociation along a streamline was calcu- 
lated by graphical integration of a plot of da/ds against s, where da/ds is 
given by the equation 

g = 2.41 X 10 20 ^)(|)(0.25 - a) (X) 

which is the first term in the expression for da/dt obtained from equation (10) 
of reference 19 with (g - a)(l + a) replaced by (O.25 - a) and with all the 




Figure 9.- Ratio of [o] to [n] as a function of [02] behind a normal shock in air. (Data 

from ref . IT- ) 

l 

multiplicative constants combined into one. The second term in the cited equa- 
tion from reference 19 accounts for three-body recombination of atoms and is 
negligible at the air densities appropriate to the problem considered here. The 
| dissociation rate kq was taken from equation (13) of reference 19- The quanti- 
ties u, p, T, and a were local streamline values. 

As was mentioned previously, the pressure variation along a streamline was 
assumed to be identical to that calculated for equilibrium inviscid flow. The u 
and T variations were found by alternately expanding the gas at constant com- 
position and dissociating it at constant pressure. This was done in small steps 
and the corresponding reaction path in the enthalpy-temperature plane was deter- 
mined through the use of the equation of state and the conservation equations for 
mass, energy, and momentum. In order to perform these operations the a, s rela- 
tion had to be specified. 

In order to obtain a relation between a and s, equation (l) was integrated 
by means of an iteration procedure in which the output a, s relation was fed back 
into the procedure as an input. The initial a, s relation was assumed to be 
a = 0 for all values of s. Fortunately, the process converged rapidly. 

The reaction rate on the body streamline was such that equilibrium was 
attained in a distance of the order of the standoff distance of the shock from 
the nose. Recombination processes were so slow that the remainder of the body 
streamline was considered to be chemically frozen. Along most of the other 
streamlines chemical equilibrium was never attained; the rate simply slowed to 
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essentially zero because of the dropping temperature and density in the expanding 
flow. Downstream of the freeze point on the streamlines the chemical composition 
was assumed to be constant. 


In order to obtain an estimate of the electron concentration at the freeze 
point, equilibrium was assumed for the H + 0 ^ N0 + + e - reaction at the local 
temperature. The rate given for this reaction in reference l8 indicates that this 
should be a good approximation. 


A calculation of the effect of error in [d] /[if] on the equilibrium electron 
concentration revealed that, if M/M should be 50 to 1 at the freeze point 
instead of 5 to 1, the electron concentration should be increased by a factor of 


( 10) 1 / 2 = 3 . 16 . Also, underestimation of the freeze-point temperature by 30°° K 
(which would be caused by neglecting the heat of formation of 10 percent NO in 
the mixture) would require that the electron concentration be increased by a fac- 
tor of about 1.15 at 5*000° K. Taken together these effects indicate that the 
electron concentration may be underestimated by a factor of about 3 - 6 . 


Finite-rate recombination of the electrons was assumed to be controlled by 
the dissociative recombination process N0 + + e" -> N + 0 and was evaluated by 
integrating the equation 

a[.-] 

dt L J 

along the streamlines. The rate k r 
taken to be equal to £e“J . 

Frozen downstream values of velocity, temperature, and density were used 
until the streamlines entered the boundary layer. Boundary- layer values were used 
for the remaining portions. 


= -k r [N0 + ] [>] (2) 

was taken from reference 18 and Fn0 + 1 was 


BOUNDARY - LAYER TREATMENT 


Boundary- layer corrections were applied only to the calculations in which 
complete equilibrium was assumed and to those in which finite-rate recombination 
of ions was assumed. Both calculations were made by using the similar solutions 
for laminar equilibrium boundary layers of reference l 6 and then matching the 
resulting profiles to the inviscid flow to account grossly for vorticity inter- 
action. Briefly, the procedure began with the calculation of velocity, tempera- 
ture, and density profiles in the boundary layer at representative axial body sta- 
tions. At each of these stations the flow properties at the edge of the boundary 
layer (y = 5) were assumed to be the same as those at the body surface (y = 0) in 
the inviscid flow field. 

The read-gas method of reference 16 was applied under the assumptions of 
equilibrium thermal properties, Lewis number of 1, and pressure gradient of 0. 
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Because of the last assumption, no stations closer to the stagnation point 

than = 0 A were used. Hence, for complete equilibrium the resulting boundary- 
hn 

layer properties should be closely determined by this method, if one excludes that 
aspect of the problem due to free- stream shear which must be treated crudely in 
any case. 

For the frozen flow the same method was used, but the inviscid surface prop- 
erties, which also are the boundary- layer edge properties, were determined from 
the frozen downstream flow field. The resulting boundary layer can be described 
as being in pseudoequilibrium, since the thermal and transport properties in ref- 
erence 1 6 are for equilibrium air. 

The flow properties used in the boundary- layer calculations at y = 0 and 
at y = 5 are given for several body stations in table II. The viscous and 
inviscid profiles of velocity, temperature, and density at a given station were 
then merged by a matching procedure which favored boundary- layer values close to 
the body and inviscid flow values far from the body. The center of the transi- 
tion region from inviscid to viscous flow solutions was arbitrarily selected to 


TABLE II.- FLOW PROPERTIES AT BODY SURFACE AMD AT EDGE OF BOUNDARY LAYER 


X 

Equilibrium flow 

Frozen 

flow 

Dn 

y - 0 

y = 5 

y = 0 

y = 5 



u, ft-sec 

-1 


0.40 

0 

10,089 

0 

8,076 

1.00 

0 

11,047 

0 

8,556 

2.00 

0 

11,474 

0 

8,868 

5.40 

0 

11,850 

0 

9,059 

11.50 

0 

11,786 

0 

9,011 


T, °K 


0.40 

400 < T < 800 

4,4iO 

400 < T < 800 

2,619 

1.00 



4,o4o 
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be that point (value of y) in the boundary-layer solution at which the velocity 
begins to drop off sharply from the value at the edge and is referred to herein 
as a matching point. In the neighborhood of this matching point the velocity, 
temperature, and density profiles were joined by a smooth fairing from the one to 
the other. Figure 10 illustrates the method and some typical profiles. 

Streamlines were followed into the boundary layer by using a mass-flow 
parameter w to identify streamlines, where w is defined by the relation 

r y 

w(y) = 2*r / pu dy 

J 0 

The concept is adopted that the mass flow interior to a surface of revolution 
generated by rotating a streamline about the body axis is the same at all axial 
locations on the body. The concept of flow streamlines as particle paths becomes 
invalid in the presence of diffusion, which should be present in a reacting bound- 
ary layer. However, in lieu of a more proper approach and recognizing that the 
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(a) Equilibrium flow, 2L = 5 . 4 . 


Figure 10.- Typical boundary- layer profiles. 


l8 




(b) Frozen downstream flow. -- = 5.i. 

D n 

Figure 10.- Concluded. 


results will be increasingly poor toward the body surface, "streamlines" were 
followed into the boundary layer by tracing the locus of constant w in the 
X / Dn > y/A plane. The streamlines found in this way are shown in figure 11. 


CONSTRUCTION OF 



AND v PROFILES 


Shock- layer profiles of electron concentration and collision frequency were 

constructed at axial stations = 0.42, 1, 2, and 11.5- Representative 

D n 

samples of these profiles are shown in figures 12 and 13 . These profiles are 
slightly different from the profiles given in references 12 and 13 for the same 
body and flight conditions. The changes represent improved knowledge of the 

H and v variations based on more detailed calculations which considered more 
streamlines than were used for the former calculations. 
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Figure 11.- Flow "streamline" paths in boundary- layer region. 


Except for the two profiles involving boundary-layer corrections, the [e - ^] 

profiles are simply cross plots of streamline calculations like those shown in 
figure 8. For EQBL the EQ profile was modified in the boundary- layer region by 
calculating the equilibrium electron concentration corresponding to the local 
temperature and density in the boundary layer. 

The reasoning for FDFEBL is a little more involved. As for the other pro- 
files the process was begun by cross-plotting points from streamline calculations. 
In the boundary-layer region the properties along streamlines were corrected for 
boundary- layer effects. However, too few streamlines were calculated to define 

the sharp dropoff in the [e“] profiles which occurs close to the body surface 

and, in contrast to EQ, [e - ] does not depend uniquely on local values of 
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Figure 12.- Electron concentration profiles. — = 

D n 


temperature and density. Therefore, the transition through this part of the 
boundary layer was determined by the following procedure: 

d[e~] r -.2 

(l) Let equation (2) be expressed by the formula — — - — = -k r [e j , where 
k r is a function of temperature only. 


ds 


(2) Divide this equation by 


u 


to obtain 



4xl0 8 



Figure 13- - Collision frequency profiles. = 5-4. 

% 


( 3 ) Choose a characteristic reaction distance based on the idea that, if the 
reaction rate were held constant, j^e“] would vanish in unit distance. In equa- 
tion form this is expressed as 

fJM = _r--i 

ds ~ L e J initial 

The definition of the characteristic reaction distance is arbitrary. Its purpose 

d[e-] 

is to establish. an order of magnitude value of which would be able to 

ds 

cause appreciable decay of the electron concentration in the flow time over the 
body. 


(4) Then, if the ratio of k r to u is held constant in equation (3)^ the 
rate at which electrons disappear decreases very quickly as ^e^j becomes smaller. 

For example, a decrease in [e“] by a factor of 10 would cause the local charac- 
teristic reaction distance to increase by a factor of 100. Thus, £e”] has a 

strong dependence on the local value of k r/ /u. Reaction equations with this 
behavior are described as " self- limit ing. " It should be noted that, except for 

d[e-] 

the dependence on k r /u, ^ depends only on the magnitude of the electron 

concentration. Thus, for comparable values of kp/u, the electron concentration 

is about the same in any flow field in which electron recombination is taking 
place according to equation (3)« 
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( 5 ) As a result, the following equations are based on the idea of maintaining 
a local characteristic reaction distance of 1 foot: 



or 



Electron concentration profiles defined by equation ( 4 ) approach zero at the 
wall, vary in a reasonable way through the boundary layer, and come close to 
joining smoothly onto profiles determined from streamline calculations. Exact 

junctions are obtained by applying constant multipliers to values of £e“J 
obtained from equation ( 4 ). The section of the FDFRBL curve in figure 12 between 
y = 0 and the point corresponding to streamline 0 was obtained in this way. 

Note that this section of the profile is independent of any reasonable choice for 
the characteristic reaction distance, since the shape is uniquely determined by 
the variation of u/k r , and the arbitrary multiplier compensates for different 
choices of the characteristic length. 

The collision frequency was evaluated from the relation 

v = i.6 x io 12 — = 4.52 x io 2 7 Be. 

i/t 

1 

t and typical results are shown in figure 13 . This expression is based on the 
assumption of a constant velocity-averaged effective collision cross section of 

Q = 3.54 x 10“ 16 cm^. The value of Q is based on the electron-molecule cross 
sections for O 2 and N 2 given in reference 20 (p. 206) and on the electron-atom 
cross section for 0 given in reference 21. (The cross sections for N and 0 
were assumed to be the same.) 

CALCULATION OF TRANSMITTED AHD REFLECTED SIGNAL STRENGTH 


The propagation problem was assumed to be equivalent to the calculation of 
the energy transmitted and reflected by a plane-parallel slab of plasma in which 

[e“J and v vary only in the propagation direction. Monochromatic plane waves 

at normal incidence impinge on one side of the slab. Reference 12 shows that this 
■ problem reduces to the solution of the differential equation 

+ k 2 E = 0 (5) 

ay 
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where E is the complex electric vector, y is distance measured in the propa- 
gation direction from the plasma boundary on the incidence side, and k is the 
complex wave number. 

A number of assumptions are implicit in this formulation of the problem. 
Among these are: (l) the plasma properties are not affected by the electromag- 
netic radiation, (2) detuning of the antenna in the presence of plasma is neg- 
lected, ( 3 ) variations of the plasma properties in directions perpendicular to 
the propagation direction are neglected, (4) curvature of the plasma is ignored, 
and ( 5 ) the boundary condition at the body surface may not be realistic. 

The plane slab model was used because it led to tractable solutions of the 
propagation equations. Because of the plane wave assumption, calculated values 
of signal attenuation are more likely to be correct for an X-band open waveguide 
antenna than for a slot antenna operating in the VHF range, such as the one on 
the RAM A1 vehicle. However, there is no restriction on the steepness of plasma 
gradients in the nonhomogeneous layer. 

The solution of equation ( 5 ) subject to appropriate boundary conditions was 
carried out by machine integration by using the Runge-Kutta procedure. Equa- 
tion ( 5 ) was reduced to the following set of simultaneous differential equations 
in real variables: 


where 


^4 + k 0 2 (fV - gW) = 0 

ay 

+ k 0 2 (fW + gv) = 0 

dy^ 


E = f + ig 


k 2 = k 0 2 (V + iW) = (a + ip) 2 




2k 



The symbols k Q , a, 0, v, cx>, 
and oop have their customary 
meanings in equations of this 
kind. The phase parameter a 
and the attenuation parameter 0 
are related to V and W as 
follows: 



Figure 14- shows a typical 
variation of V and W through 
a plasma layer. Continuity of 
V and W at the boundaries is 
not required and will not be 

present if [e - ] is discon- 
tinuous at the boundaries. In 
this figure the incident wave is 
indicated by the unit amplitude 

plane wave e^°^. The reflected 
and transmitted waves are also 
indicated as plane waves with the 



Figure l 1 )-.- Variation of V and W in the shock 

layer. FDFRBL: — = 5.1).. 

D n 


complex amplitudes Er and Erp, respectively. 


The values of V and ¥ in mediums 0 and 2 are taken to be equal to the values 
in vacuum. The waves are assumed to be simple harmonic in time with frequency co. 


Continuity of E and E' at the boundaries requires that the following 
relations be satisfied: 


E 0a = 1 + I Er| e k a = E 


'R 


E, 


i0a = 

3a = ik o ~ iko^e 1 ^ = E-[ 

E ib = | E t Je^(k° A+ ^tJ 

E k - ik 0 |ii r |. 1 ( k ° a+ i iS b) , 

The propagation equations cannot be numerically integrated in the propaga- 
tion direction because Ep is unknown. The integration can be performed in the 


E la 

(6) 


(7) 

11 

c? 

(8) 

)_E' 

“ Ji 2b 

(9) 
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opposite direction, however, 

. , ifkoA+^h) 

divided by Ej e ' , 


if the equations for the boundary conditions are 
as follows: 


%a - 1 + %f° ^ a - E la 
®0a ' = Eia 


E lb ~ 1 “ E 


2b 


E lb = ik o = E 2b 


where 


I = 


e r - 


1 %!' 


l&T ! 


Eji 


The quantities needed for the calculation of the transmitted and reflected 
signal strength are | Ej | and | Er | from equations (6) to (9). These are related 

to 1 1 j and J Er | as follows : 

W-rfi 


E R | = 


Kl 


The integration was begun at y = A where f = 1 , f = 0 , g = 0 , g = k Q , 
and was continued to y = 0 to find the values of f, f g, and g' at 
y = 0 . These values were then substituted in the following formulas to find jl 
and I Er I : 


n.\ 2 


,.\2 


9 'iH( k ° J + r ) + ( k °® - 


%| - irll( ko? - g ') 2 + (*o 5 + ?') 


-,.l2 
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The final results on reflected and transmitted signal strength were obtained 
in terms of signal power loss in decibels. The transmission loss due to reflec- 
tion from the plasma is 

(SL) h = 10 log ^1 - |e r | 2 ) 

The transmission loss due to both reflection and attenuation is 

(SL) t = 20 log | Ep | 

The calculated transmission loss for various flow assumptions is given in 
figure 15 . The ratio of reflection loss to loss caused by both reflection and 
attenuation is given in figure 16 . 

The normal incidence program described has been generalized to cover plane 
waves incident at any angle on a nonhomogeneous slab . The new program was 
reported in reference 12 and was applied to a study of plasma reciprocity (inter- 
change of incidence and transmission sides of the plasma) in reference 22 . 




Figure 15 • - Signal loss as a function of 
axial distance from the nose. 


Figure l6.- Fraction of signal lost by- 
reflection . 
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CONCLUDING REMARKS 


Even on the basis of the crude models used in this study, the indication is 
clear that adequate description of plasma properties for the purpose of calcu- 
lating radio signal attenuation will in general require that allowance be made 
for (l) finite rate chemistry and ionization and ( 2 ) the presence of a boundary 
layer. 

The plane slab model used for calculation of signal attenuation is more 
likely to give good results when the signal wavelength is small compared to body 
dimensions, which was not the case for the VHF slot antenna on RAM Al. Thus, the 
numerical values of signal attenuation given herein are questionable and should 
be recalculated as soon as a program is available which is based on a geometry 
more closely approximating the actual case. Other areas in which improvements 
might be made are consideration of antenna-plasma interactions and proper defini- 
tion of boundary conditions at the boundary near the antenna. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 10 , 1963. 
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